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We introduce a Gaussian version of the entanglement of formation adapted to bipartite Gaussian 
states by considering decompositions into pure Gaussian states only. We show that this quantity is an 
entanglement monotone under Gaussian operations and provide a simplified computation for states 
of arbitrary many modes. For the case of one mode per site the remaining variational problem can 
be solved analytically. If the considered state is in addition symmetric with respect to interchanging 
the two modes, we prove additivity of the considered entanglement measure. Moreover, in this case 
and considering only a single copy, our entanglement measure coincides with the true entanglement 
of formation. 

PACS numbers: 03.67.Mn, 03.65.Ud, 03.67.-a 



I. INTRODUCTION 

One of the main novelties of Quantum Information 
Theory is to consider entanglement no longer merely as 
an apparent paradoxical feature of correlated quantum 
systems, but rather as a resource for quantum informa- 
tion processing purposes. This new point of view natu- 
rally raises the question regarding the quantification of 
this resource. How much entanglement is contained in 
a given state? For pure bipartite states there is, under 
reasonable assumptions, a simple and unique answer to 
this question, namely the von Neumann entropy of the 
reduced state [l|, II, 13] • F° r mixed states there are several 
entanglement measures which can be distinguished 
due to their operational meaning and mathematical prop- 
erties. Such a measure should be non-increasing under 
mixing as well as under local operations and classical 
communication (LOCC), and it should return the right 
value for pure states. The largest measure fulfilling these 
requirements is the Entanglement of Formation Ep |5( • 
Operationally, it quantifies the minimal amount of entan- 
glement, which is needed in order to prepare the state by 
mixing pure entangled states. It is therefore defined as 
an infimum 

E F {p) = mi\Y jPk E{^ k )\ p = Y,Pk\^k)^k\\ (1) 

I k k ) 

over all (possibly continuous) convex decompositions of 
the state into pure states with respective entanglement 
E(V) = 5(tr B [|*)(*|]), where S(X) = -trpflogX] is 
the von Neumann entropy. By its definition calculating 
Ep is a highly non-trivial optimization problem, which 
becomes numerically intractable very rapidly if we in- 
crease the dimensions of the Hilbert spaces. Remarkably, 
there exist analytical expressions for two-qubit systems 
|(| as well as for highly symmetric states Q . 

Recently, Ep was calculated for the first time for con- 
tinuous variable states, namely for symmetric Gaussian 
states of two modes Q . In general, Gaussian states are 
distinguished among other continuous variable states due 



to several reasons. Experimentally, they are relatively 
easy to create and arise naturally as states of the light 
field of a laser (cf.[lfj|) or in atomic ensembles interact- 
ing with light [ll| . For this and other reasons they play a 
more and more important role in Quantum Information 
Theory (cf. [13). 

Theoretically, despite the underlying infinite dimen- 
sional Hilbert space, they are completely characterized by 
finitely many quantities - the first and second moments 
of canonical operators. Moreover, they stand out due to 
several extremal properties |27j. In fact, the calculation 
of Ep for symmetric two-mode Gaussian states depends 
crucially on the fact that for given "EPR-correlations" 
two-mode squeezed Gaussian states are the cheapest re- 
garding entanglement. This implies that in this particu- 
lar case there is a decomposition in terms of pure Gaus- 
sian states, which is optimal for Ep in (|T]l. 

On the one hand, this raises the question, whether 
this is generally true for all Gaussian states. On the 
other hand one may, motivated by the operational in- 
terpretation of Ep, restrict Eq. JQ) to decompositions 
into Gaussian states from the very beginning. After all, 
Gaussian states arise naturally, whereas the experimen- 
tal difficulties of preparing an arbitrary pure continuous 
variable state are by no means simply characterized by 
the amount of its entanglement. For these reasons we will 
in the following investigate the Gaussian Entanglement 
of Formation Eq to quantify the entanglement of bipar- 
tite Gaussian states by taking the infimum in Eq. 
only over decompositions into pure Gaussian states. 

This article is organized as follows: in Sec. [H] we recall 
basic notions concerning Gaussian states. Sec. IIIII de- 
fines the Gaussian Entanglement of Formation and pro- 
vides a major simplification concerning its evaluation for 
bipartite Gaussian states of arbitrary many modes. In 
Sec. IIVI we prove that Eq is indeed a (Gaussian) entan- 
glement monotone, in the sense that it is non-increasing 
under Gaussian local operations and classical communi- 
cation (GLOCC). The case of general two-mode Gaussian 
states is solved analytically in Sec. El The special case 
of symmetric Gaussian states, for which it was proven 
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m @ that E G = E F is discussed in detail in Sec. IVII 
where we give an alternative calculation of Eq, which is 
in turn utilized in Sec. IVIII in order to prove additivity 
of Eq for this particular case. Finally, Sec. IVIlH applies 
the measure to some examples which arise when a two- 
mode squeezed state is sent through optical fibers. The 
appendix proves a Lemma about decompositions of clas- 
sical Gaussian probability distributions. 



II. GAUSSIAN STATES 



Consider a bosonic system of n modes, where each 
mode is characterized by a pair Qk,Pk of canonical 
(position and momentum) operators. If we set R = 
{Qi, Pi, ... , Qm Pn) the canonical commutation relations 
are governed by the symplectic matrix 



fc=i 



1 

-1 



(2) 



via [Rk,Ri] = i<Tki- A state is called a Gaussian state if 
it is completely characterized by the first and second mo- 
ments of the canonical operators Rk in the sense that the 
corresponding Wigner function is a Gaussian. Utilizing 
Weyl displacement operators W% = e 1 ^ , the first mo- 
ments dk = tr [pRk] can be changed arbitrarily by local 
unitaries. Hence, all the information about the entangle- 
ment of the state is contained in the covariance matrix 
(CM) 



jkl = tr [p{R k - d k l, Ri - 



(3) 



where {•, •}+ denotes the anti-commutator. By definition 
the matrix 7 is real and symmetric, and due to Hcisen- 
berg's uncertainty relation it has to satisfy 7 > ia. For 
pure Gaussian states we have det(7) = 1 or, equivalcntly, 
(a 7 ) 2 = -l. 

In the following we denote the density operator corre- 
sponding to the Gaussian state with covariance matrix 
7 and displacement vector d by P( 7 ,d). If the latter is a 
bipartite state, its tensor product structure corresponds 
to a partition of the n modes into two subsets. 

An important decomposition of P( 7 ,o) m t° pure Gaus- 
sian states is given by 



7*.d-£) e 



(4) 



where j p < 7 is the covariance matrix of a pure Gaussian 
state. Since displacements of the form P( 7 .o) l ~ > PH,d) are 
local operations, Eq. (@J tells us that starting with P( 7 ',o) 
we can obtain every Gaussian state with CM 7 > 7' by 
means of LOCC operations. 



III. GAUSSIAN ENTANGLEMENT OF 
FORMATION 

We define the Gaussian Entanglement of Formation 
Eq for a bipartite Gaussian state P( 7) d) by 

E G (p(~f,d)) = inn / p{d-y p ,dD) E(p {jp , D) ) (5) 



P(n,d)= J P-( d lp,dD)p {lp . D) j, (6) 

where the infimum is taken over all probability measures 
p characterizing convex decompositions of pij,d) m t° pure 
Gaussian states p( 7 u), and E{p( rt p ,d)) is the von Neu- 
mann entropy of the reduced state. For pure n x n mode 
Gaussian states this quantity can be readily expressed in 
terms of the symplectic eigenvalues of the reduced CM. 
Denote these eigenvalues by a^,k — 1, . . . , n. We have 
aj. > 1 and define > by = coshrv Then 



(7) 



where 



H(r) — cosh 2 (r) log 2 (cosh 2 r) — sinh 2 (r) log 2 (sinh 2 r). 

(8) 

To obtain this expression note that any pure nxn Gaus- 
sian state is locally equivalent (via unitary GLOCC) to 
the tensor product of n two-mode squeezed states with 
squeezing parameters r k 01- For each tensor factor the 
entanglement is given by the above formula. Since the 
symplectic spectrum of the reduced CM is invariant un- 
der local unitary GLOCC the r k can be computed di- 
rectly from 7 P as described above. 

The integrals in Eqs. 10 EJ ar e taken over the space 
R n of displacements and over the set of admissible pure 
state covariance matrices. The following proposition tells 
us that it is sufficient to only consider measures p, which 
vanish for all but one covariance matrix: 

Prop. 1 The Gaussian entanglement of formation for a 
bipartite Gaussian state P( 7 ,d) is given by 



E ( 



(P( 7 ,d)) = inf \E(p hpfi) ) I 7 P < 7}, 



(9) 



where the infimum is taken over pure Gaussian states 
with CM 7 P . 

Proof: The proof can be divided into three steps: 

(i) The problem can be reformulated in terms of classi- 
cal Gaussian distributions, by considering Wigner func- 
tions instead of density operators. This is formally 
achieved by taking the trace of the decomposition in 
Eq. © with t he p hase space displaced parity operator 
T>t = WfVWl 0. Then 

tr [PiPh-d)] = bi~* ex P [ - \(<r£ + dfj- 1 ^ + d) 

(10) 
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is up to a normalization factor equal to the Wigner func- 
tion of P( 7!< j), which in turn completely determines the 
state. 

(ii) All the states P( 1p ,d) contributing to Eqs. (0 
must have smaller covariance matrices 7 P < 7, i.e. 



{{1p,D) 



{x, (7 - 1p)x) < 



}) 



(11) 



The idea of the proof is that the tails of a Gaussian with 
CM that is too large would give rise to an increasing and 
in the end overflowing contribution if we only move far 
enough away from the center. This is mathematically for- 
malized in Lemma |21 stated and proven in the appendix. 
To apply Lemma|31via Eq. (jlOJI to Eq. © we need in ad- 
dition, that the inverse is operator monotone on positive 
matrices 7 > 0. 

(iii) Assume that p is a measure corresponding to an 
optimal decomposition of p^^d) giving rise to the infimum 
in Eq. ©. Then 



E G (P( 7 , d )) = / Kdlp,dD)E(p hptD) ) (12) 

■>7p<1 

> inf{£( P(> , 0) ) I 7P <7}- (13) 

However, by using a Gaussian decomposition of the form 
in Eq. (0J we know that equality in Eq. (|13fl can be 
achieved for a measure p which is Gaussian in D and 
a delta function with respect to 7 P . □ 

Proposition^considerably simplifies the calculation of 
Eg, since the optimization is reduced from the set of all 
possible decompositions to the set of pure states satisfy- 
ing the matrix inequality 7 p < 7. Before we proceed to 
calculate Eg analytically for the two-mode case, we will 
show that Eg is indeed a proper entanglement monotone. 
To confine the argument to CMs (rather than density ma- 
trices) we make use of Proposition ^ 



IV. MONOTONICITY UNDER GAUSSIAN 
OPERATIONS 

For Eg to serve as a good Gaussian entanglement mea- 
sure it should not increase under GLOCC. That this 
is the case is quickly seen u sing the characterization of 
Gaussian operations given in |l5lll6| . There it was shown 
that the change of the CM 7 of a Gaussian state under 
Gaussian operations takes the form of a Schur comple- 
ment: 



G( 7 ) =T 1 -f 1 



1 



r 2 + 7 



fJ 2 = S(t + 0© 7 ). (14) 



Here 



Ti r 12 
rf 2 r 2 



is the 4n x 4n CM of the state characterizing the opera- 
tion G and T > is the CM of the partially transposed 
state 28]. 

For positive matrices A > B implies S(A) > S(B) (cf. 
[T7L p.472]). Consequently, if 7 > j p , then the trans- 
formed CM fulfills G( 7 ) > G( 7p ). 

Every GLOCC can be decomposed into a pure GLOCC 
G p mapping pure states onto pure states, and the addi- 
tion of classical Gaussian noise. This decomposition can 
easily be shown using the above mentioned ordering of 
the Schur complements for ordered matrices. The de- 
composition then reads G(7 P ) = G P (7 P ) + P, where the 
noise is characterized by some positive matrix P > 0, 
which is usually state-dependent. Therefore we have 
G(7) > G(7 P ) > G P (7 P ) and since the latter CM cor- 
responds to a pure state, which can be obtained from 7 P 
by a local Gaussian operation, its entanglement is cer- 
tainly smaller than that of 7 P 0|- It follows that Eg 
cannot increase under GLOCC. 



V. THE GENERAL TWO MODE CASE 

Now that we have assured that Eg is a good measure of 
entanglement in a Gaussian setting, we set to compute it 
for the case of two Gaussian modes in an arbitrary mixed 
state. The CM 7 of any two-mode Gaussian state can be 
brought to the normal form [T^. I20I ] 



( n a 






k p \ 


u 









C„ © G v 



(15) 



with k q > |A: p | by local unitary Gaussian operations. 
The block structure corresponds to a direct sum of po- 
sition and momentum space, i.e., we have reordered 
R = (Qii Q2, Pi, P 2 ). Since the normal form in Ea. lflTf is 
unique the parameters (n a , rn,, k q , k p ) provide a complete 
set of local invariants. 

The first step towards calculating Eg for these states 
is to show that there is always a pure state 7 p , which is 
optimal for Eq. © and has the same block structure as 
7 m Eq. (P3). To this end we will first provide a general 
parameterization for pure state CMs, which accounts for 
the direct sum with respect to configuration and momen- 
tum space: 

Lemma 1 A real symmetric matrix 7 P is the covariance 
matrix of a pure Gaussian state of n modes iff there exist 
real symmetric nxn matrices X and Y with X > such 
that 



7 P = 



X XY 
YX YXY + X- 1 



(16) 



where the block structure corresponds to a direct sum of 
configuration and momentum space. 

Proof: A covariance matrix 7 P corresponds to a pure 
Gaussian state iff (7 P cr) 2 = 1. If we write 



7 P 



X C 
C T D 
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FIG. 1: For any two-mode Gaussian state with CM 7 = C q ffi 
C p the Gaussian Entanglement of Formation is given by the 
entanglement of the least entangled pure state with CM 7 P = 
X X" 1 which is such that C" 1 < X < C q . Moreover, the 
optimal X can be shown to lie on the rim of the intersection of 
the forward and backward cones of CZ 1 and C q respectively. 



with X, D > 0, then this is equivalent to 

XD = 1 + C 2 , 
DC = (DC) T , 
CX = (CX) T . 



(17) 
(18) 
(19) 



Eq. (H) implies that Y := X~ X C = (X- 1 C) T = Y T 
is indeed symmetric. Moreover, Eq. I|17|) leads to D = 
X~ x (l + C 2 ) = X^ 1 + YXY. Hence, every covariance 
matrix of a pure Gaussian state is of the form in Eq. (|16fl . 

Conversely, every such matrix j p with X > is posi- 
tive definite and has symplectic eigenvalues equal to one 
since the spectrum of — (7 p er) 2 is the symplectic spec- 
trum squared of 7 p . Thus every matrix j p is an admissi- 
ble covariance matrix corresponding to a pure Gaussian 
state. □ 

The covariance matrix in the normal form of Eq. i|15|) 
only contains terms which are quadratic in the mo- 
menta but has no linear contributions. This implies that 
the state remains invariant under momentum reversal 
P ^> —P and since this can be interpreted as complex 
conjugation, the respective density operator is real (in 
position representation) . 

Eq. Ijl6|l gives the covariance matrix of a pure Gaussian 
state with respective wave function 

#0) = \irX\-i cxp [ - ^(X- 1 - iY)x] , (20) 

which in turn becomes real if Y = 0. The following 
Lemma shows that for two- mode states we can, in fact, 
restrict to these real pure states in the calculation of 
B G (7): 

Lemma 2 Let 7 = C q © C p be the covariance matrix of a 
two-mode Gaussian state. Then there exists a pure state 
with covariance matrix of the same block structure which 
minimizes Eq for 7. 



Proof: We will show that for every j p (X,Y) < 7 of 
the form in Eq. (|16|) the covariance matrix j p (X, 0) leads 
to an improvement for Eq. First note that the block 
structure of 7 implies that 7 > 7 P (X, ±Y") and thus 



1 



7 > ~[ lp (X,Y) 



lp{X,-Y)] 
= lp (X 7 0) + 0®(YXY) 



(21) 
(22) 
(23) 



Therefore ^p{X, 0) is an admissible covariance matrix for 
the Eg optimization problem. 

In order to show that j p (X, 0) is less entangled than 
7 P (X, Y) we make explicit use of the assumption that we 
deal with two-mode states. In this case the entanglement 
is a monotonous function of the determinant of the re- 
duced covariance matrix. The difference of the respective 
determinants can be calculated straight forward and it is 
given by 



det [i A) (X,Y] 
= Y? 2 det [X] 



- det [7 
> 0, 



(24) 
(25) 



which completes the proof. □ 

According to Lemma El the remaining task for calcu- 
lating Eq is to find the CM j p = X © X~ x which has 
minimal entanglement under the constraint that 



c- 1 <x<c q . 



(26) 



This inequality has a simple graphical depiction stem- 
ming from the fact that the set of positive semi-definite 
matrices X satisfying an inequality as, e.g. C q — X > 0, 
form a cone, which is equivalent to the (backward) light 
cone of C q in Minkowski space: if we expand a Hermitian 
2x2 matrix in terms of Pauli matrices (and the identity) , 
the expansion coefficients play the role of the space-time 
coordinates and the Minkowski norm is simply given by 
the determinant of the matrix. Hence, by Eq. l|26|l X has 
to lie in the backward cone of C q and in the forward cone 
of C" 1 (see FigCJ). 

Instead of minimizing the entropy of the reduced state 
under this constraint, we may as well minimize the de- 
terminant of one of the local covariance matrices 



m(X) 



1 



dctX 

XniX- 1 ), 



1 



(x- 



detpsT- 1 ) 



(27) 
(28) 



since, as already stated, this is a monotonously increasing 
function of the entanglement. Thus we have to find 



min {m(X) \C p 1 <X <C q ) 



(29) 



over the real, symmetric 2x2 matrices X. 

In fact, for the optimal X both inequalities have to be 
saturated, i. e. 



det (C q -X) = det (X - C^ 1 ) = 0. 



(30) 
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In order to see this assume we are given a matrix X with 
Cp 1 < X < C q . Then we can decrease the value of 

m(X) with a matrix X := X + el by increasing e > 
until C q — X is of rank 1. However, by Eq. 127|) the same 
argument holds for X^ 1 < C p . 

To depict it geometrically again, the optimal X has to 
lie on the rim given by the intersection of the backward 
and forward cones of C q and C" 1 respectively. Hence, 
we have reduced the number of free parameters in the 
calculation of -Eg (7) to one angle, which parameterizes 
the ellipse of this intersection. 

For every explicitly given CM 7 minimizing m(X) on 
this ellipse is now straight forward. Writing down the 
resulting value for Eq in terms of the general parame- 
ters {n a ,nb,kq,kp) of Eq. (|15fl leads, however, to quite 
cumbersome formulae involving the roots of a forth or- 
der polynomial. Since not much insight is coming out 
of these expressions we refrain from writing them down 
explicitly and continue with discussing the special case 
n a = rib f° r which we obtain a simple formula for Eq. 

Nevertheless, for an arbitrary but explicitly given 7 
the remaining variation under the constraint in Ea. l|3UI) 
is a simple exercise which can be solved analytically with 
the help of any computer algebra program. For some 
examples, see Sec. IVIIII 

VI. SYMMETRIC STATES 

Symmetric two-mode Gaussian states with CM of the 
form in Eq. (|15|l with n a = nb = n arise naturally when 
the two beams of a two-mode squeezed state are sent 
through identical lossy fibers [5l| (see also Sec lVlll|) . The 
entanglement of formation Ep of these states was calcu- 
lated in and it was proven that a decomposition into 
Gaussian states gives rise to the optimal value. Together 
with the obvious fact that Eq is an upper bound for Ep 
this implies that Eq = Ep in this case. Since the calcu- 
lation of Ep is however quite technical and in order to 
make the present article more self-contained, we provide 
in the following a simpler way to obtain Eq. 

In principle we could utilize the general results of the 
previous section, which simplify greatly for the symmet- 
ric case. However, we give an alternative proof and re- 
duce the result to the fact that the optimal j p in Eq. © 
has the same logarithmic negativity [221 ] as 7. A similar 
argument is used in Sec. IVIll to prove additivity of Eq. 

Prop. 2 (Eq for symmetric states) For symmetric 
lxl Gaussian states, i.e. states whose CM 7 is char- 
acterized by local invariants (n,n, k q , k p ), the Gaussian 
entanglement of formation is given by 

E G (rf) = H(r ), (31) 

where the minimum two-mode squeezing required is given 
by 

r = ^\n[{n - k q ){n + k p )} (32) 



and H(r) is defined in Eq. Q). 

Proof: First, instead of 7 we consider the locally equiv- 
alent CM 7' which is obtained from 7 by squeezing |2^| 
both Qa and Qb by A = [(n + k p )/(n — kq)] 1 / 4 . Clearly 
the CM 7' has the same Eq as 7. It is straightforward 
to check that the pure two-mode squeezed state with 
two-mode squeezing parameter ro and corresponding CM 
7(j"o) is indeed smaller than 7'. 

That there can be no pure state j p with less entan- 
glement satisfying j p < 7' follows from the monotonic 
dependence of pure state entanglement on the two-mode 
squeezing parameter: any pure two-mode Gaussian state 
is locally equivalent to a two- mode squeezed state 7(7*) 
and its entanglement is given by H(r). An important 
entanglement-related characteristic of these CMs are the 
symplectic eigenvalues of the partially transposed CM 
[22|, in particular those smaller than one. They are 
invariant under local unitary Gaussian operations and 
for the two-mode squeezed s tate given by e ± r . For the 
symmetric CM 7 they are y/(n=p k q )(n ± k p ). Thus the 
smallest symplectic eigenvalues of 7 and 7(ro) coincide. 

For positive matrices A > B implies > bk, where Ofc 
(bk) denote the ordered symplectic eigenvalues of A (B) 
jl8j . Since the ordering A > B is preserved under partial 
transposition, all pure states with less entanglement than 
7(ro) cannot possibly satisfy 7 > j p , hence ro is optimal. 

□ 

Thus for symmetric states 7 the optimal pure state 7 P 
is characterized by the fact that the smallest symplectic 
eigenvalues £1(7) and si(7 p ) of the two partially trans- 
posed CMs are identical. According to [221 this implies 
that the logarithmic negativity, Sat (7) = — ^ln[si(7)] °f 
both states is the same, i.e., in the optimal decomposition 
pure Gaussian states are mixed such that "no negativity 
is lost" in the mixing process. For non-symmetric states 
this is no longer possible and Si("/) is strictly larger than 
si(7p), i.e. more entanglement is needed to form 7 than 
required by its negativity. 

VII. ADDITIVITY 

One longstanding question about the entanglement of 
formation is if it is additive, that is whether Ep(p\®p2) = 
Ef(pi) + Ep(p2) or whether one may get an "entan- 
glement discount" when generating several states at a 
time. Here we show that for symmetric Gaussian states 
the Gaussian entanglement of formation Eq is additive. 
Since for these states Eq was shown [j| to equal Ep this 
may hint at additivity of even the latter quantity for 
Gaussian states. 

Prop. 3 (Eq is additive for symmetric states) Let 

= 1, . . . , N describe symmetric Gaussian states with 
local invariants (ni,ni,k q j,k p j) and let T — ®fL 1 'ji 
describe the tensor product of these states, then 

E g (T) = Y. E g(^)- ( 33 ) 
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Proof: Let the logarithmic negativity of the fcth state 
be r"fe, and assume r^+i < 77.. To show additivity, we use 
again that A > B > implies > bk for the ordered 
symplectic eigenvalues of A{B). 

Let T p < T be a pure N x TV-mode CM. Consider the 
partially transposed CM T. We have 



Eg(p) 



r>r p , 

which implies > s k , where {sk,k = 1,...,N} denote 
the (descendingly ordered) symplectic eigenvalues of f, 
and s' k the same for r p . 

All pure bipartite Gaussian states are locally equiv- 
alent to a tensor product of two-mode squeezed states 
|13| with (ordered) two- mode squeezing parameters tfe. 
For the following we only need to look at the ./V smallest 
symplectic eigenvalues. For these we have 



for T r 



0.5 



l/lt 



0.1 



0.2 



0.3 



0. 



0.5 



FIG. 2: Gaussian Entanglement of Formation (in units of 
ebits) for a TMSS with r = 1 after transmission through a 
lossy optical fiber at temperature r — 0. The plot shows 
Eg(p) versus transmission length I/Ia for the symmetric and 
asymmetric setting (dotted and solid line, respectively). 



and 



Sk 



for r. 



Hence T > T p implies t& > rfe, i.e., the optimal joint 
decomposition is the tensor product of the optimal de- 
compositions for the individual copies. Thus Eg is addi- 
tive. □ 

For the non-symmetric case the optimal individual de- 
composition does no longer allow tk = ffe, i.e. more en- 
tanglement than required by the logarithmic negativity 
must be expended to produce p. Therefore, the previous 
argument does not hold and the question of additivity of 
Eg remains open for general lxl Gaussian states. 



length, one or the other setting yields more entangle- 
ment for the distributed state from a given squeezing of 
the initial state. 

According to [23( the Gaussian state after transmission 
through the fiber has a CM 
















(35) 



with c'i 
and s = 



s' 
(2N th + 1) (1 



for i = 1,2 



VIII. EXAMPLES 

In this section we will apply the Gaussian Entangle- 
ment of Formation to a simple practical example. Con- 
sider a two-mode squeezed state (TMSS) with CM 

/ c -s 0\ 

c s I 

7 ~ s c ' (34) 

\ s c/ 

c = cosh2r, s = sinh2r, 

which is to be distributed between two parties by means 
of a lossy optical fiber. There are two extremal settings 
for the transmission of the state which may lead to differ- 
ent values for the distributed entanglement: The source 
could be placed either at one party's site ("asymmetric 
setting") or halfway between both parties ("symmetric 
setting"). In the former case, one mode is transmitted 
through the whole length of the fiber while the other one 
is retained unaffected, in the latter both modes are trans- 
mitted through half the length of the fiber each. It turns 
out that depending on thermal noise and transmission 



Here Ti, Ti are the transmission coefficients of the fiber 
for the respective modes and the fiber is in a thermal 
state with mean number of photons Nth, which is in turn 
related to a "temperature" r by A^h = (exp(l/r) — 1) _1 . 
Note that in quantum optical settings using optical fre- 
quencies we have r = 0. 

Depending on the setting, the transmission coefficients 
take on the values T\ — e~ l / lA ,T2 = 1 (asymmetric set- 
ting) or T\ — T2 — e~ z / ( 2 ^a) (symmetric setting) where 
I /I A denotes the total length I of the fiber in units of the 
absorption length I a- As an example, we compare the 
two settings for a TMSS with r = 1. While for tempera- 
ture r = (Fig. 0) the asymmetric setting always yields 
a higher entanglement for the final state, at finite tem- 
perature t = 1 (Fig. 21 the symmetric setting is to be 
preferred for longer ranges. 

Fig. 21 shows Eq as a function of the initial squeezing 
r and the transmission length I for the symmetric setting 
at zero temperature. As already indicated in 0, |2j| in- 
creasing the squeezing over a certain threshold has only a 
negligible effect on the transmitted entanglement already 
after a small fraction of the absorption length. 
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FIG. 3: Gaussian Entanglement of Formation (in units of 
ebits) for a TMSS with r = 1 after transmission through a 
lossy optical fiber at temperature r = 1. The plot shows 
Eg(p) versus transmission length I /I a for the symmetric and 
asymmetric setting (dotted and solid line, respectively). 




FIG. 4: Gaussian Entanglement of Formation (in units of 
ebits) for a TMSS with initial squeezing r after symmetric 
transmission through a lossy optical fiber at zero tempera- 
ture. I/Ia is the transmission length in units of the absorption 
length. 



IX. CONCLUSION 

We introduced a Gaussian version of the Entanglement 
of Formation by taking into account only decompositions 
into pure Gaussian states. On Gaussian states this is 
a proper entanglement measure in the sense that it is 
non-increasing under GLOCC operations. Moreover, it 
is an upper bound for the full Entanglement of Forma- 
tion which is tight and additive at least for symmetric 
two-mode Gaussian states. However, it remains an open 
question whether this is true for all Gaussian states. 

We have shown how to analytically calculate Eq for all 
two-mode Gaussian states and given a simple formula for 
the symmetric case. For multimode bipartite states, the 



(numerical) computation of Eq becomes rapidly more 
difficult: Although Proposition 1 still allows to replace 
the minimization over all possible decomposition into 
pure Gaussian states by the minimization over the (for 
the n x n case) 2n(2n + 1) parameter set of pure 2n-mode 
covariance matrices the further simplifications used in the 
lxl case are no longer available and the minimization 
is no longer easy. To compute Eq for interesting multi- 
mode states such as two copies on non-symmetric lxl 
states or bound entangled states j2(| new methods need 
to be developed. 
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QUPRODIS and the German National Academic Foun- 
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APPENDIX A: PROOF OF LEMMA 3 

The following Lemma about decompositions of clas- 
sical multivariate Gaussian distributions is used in the 
proof of Prop. ^ 

Lemma 3 Let 

G(A,a,x) := |A|*(27r) _ * e -^ x - a ^ A< - X ~ a \ x,ae M n 

(Al) 

be a Gaussian probability distribution with symmetric 
A > and consider an arbitrary convex decomposition of 
G{A, 0, x) into other Gaussian distributions of this form: 

G(A,0,x) = J fj,(dB,db)G(B,b,x). (A2) 

Then all the distributions contributing to this decomposi- 
tion have to satisfy B > A in the sense that 

fi({(B,b) | 3 X : x T (A - B)x > 0}) = 0. (A3) 

Proof: Let us first define a set 

S(x,e) := {{B,b) \ x T Ax- e\\x\\ 2 >x T B x y (A4) 

Integrating Eq. (|A2|I only over S(x, e) leads then to a 
lower bound on G(A, 0,x): 

G(A,0,x)> [ n(dB,db)G{B,b,x). (A5) 

JS(x,e) 

Inserting G(A, 0,x) and G(B,b,x) we can symmetrize 
inequality l|A5(l with respect to x i— > —x, which leads to 

\A\i e -? xTAx > [ n(dB,db) \B\i cosh (x T Bb) 

JS(x,e) 

e -±(x T Bx+b T Bb) _ / A gj 
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Utilizing cosh > 1 and the defining property of the set where we have of course used that Q" is dense in 1™. □ 
S(x, e) we have 

l^e-ill*!! 2 > f n(dB,db) \B\h-? bTm . (A7) 

JS(x.e) 

Note that the r.h.s. of Eq. (|AT|) does no longer depend 
on the norm of x but merely on its angular components. 
Taking the limit ||x|| — > oo implies then that S(x,e) is of 
measure zero for every e > 0. Moreover, every countable 
union of such sets is of measure zero. In particular: 

U U = (A8) 

xeQ* 1 m£N 

= {(B,b)\3x e Q",3m £ N : 

x T Ax > x T Bx + — \\x\\ 2 ) (A9) 
m ) 

= ((B,b)\3x S Q™ : x T Ar > x T Bx\ (A10) 

= {(fl,&)|3a: € R" : x T (A-B)x > o}, (All) 
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